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Abstract 

We consider a classical test particle subject to electromagnetic and 
gravitational fields, described by a Lagrangian depending on the ac- 
celeration and on a fundamental length. We associate to the par- 
ticle a moving local reference frame and we study its trajectory in 
the principal fibre bundle of all the Lorentz frames. We discuss in 
this framework the general form of the Lagrange equations and the 
connection between symmetries and conservation laws (Noether the- 
orem). We apply these results to a model, already discussed by other 
authors, which implies an upper bound to the proper acceleration and 
to another new model in which a similar quantity, called "pseudo- 
acceleration" , is bounded. With some simple choices of the fields, we 
illustrate some other interesting properties of the models and we show 
that unwanted features may appear, as unstable run-away solutions 
and unphysical values of the energy-momentum or of the velocity. 
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I Introduction. 



At present, the only well established long range relativistic fields, that allow 
a classical (non quantum) treatment, are the Maxwell electromagnetic field 
and the Einstein gravitational field. However, the existence of new long range 
fields, or of modifications of the above mentioned ones is often suggested on 
the basis of theoretical speculations or of not yet well confirmed experimental 
results. Many crucial problems of present day physics and cosmology can be 
formulated in terms of long range fields, for instance a possible failure of the 
Lorentz symmetry can be attributed to a long range nonscalar field. 

A long range field can be defined in terms of well chosen test particles. 
For instance, charged particles are needed to define the electromagnetic field 
and only spinning particles are affected by a torsion field [1] . We see that the 
investigations of long range classical fields and of test particles are intimately 
related. The aim of the present article is to give a contribution to this problem 
by discussing new kinds of test particles. 

By definition, in a test particle model one disregards the influence of the 
particle on the field, namely several important effects, as the radiated energy, 
the radiation damping, the electromagnetic contribution to the mass, etc. 
However, these models are much simpler and permit a better undertanding 
of some general features of the theory, in particular its symmetry properties. 

The test particle models may also suggest a way to introduce in the 
theory a fundamental length A of the order of the Planck length with the 
aim of taking into account approximately some effect of quantum gravity. 
For instance, as it has been proposed by several authors [2-10], an effect of 
quantum gravity should be an upper bound of the order of A -1 to the proper 
acceleration of particles (we use the convention c = 1). Note that the gravi- 
tational acceleration vanishes in a local inertial frame and a discussion of the 
maximal acceleration principle requires the presence of an electromagnetic 
field. More references are given in the report [11], dedicated to a discussion 
and a comparison of various geometric descriptions of the maximal accelera- 
tion principle. Anther effect of the fundamental length is the appearance of 
new degrees of freedom and of excited states of the particle. 

In order to study the motion of a test particle, it is useful to associate 
to every point of its world line, parametrized by the variable r, a moving 
local Lorentz frame (a tetrad) s(r) and to consider the motion of the point 
s(r) in the principal bundle S of all the local Lorentz frames. This point 
of view has been adopted in refs. [12, 13] in the framework of presymplectic 
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dynamics [14] and in refs. [15-17] starting from the balance equations of a 
field theory and following the ideas of ref. [18]. 

The method of moving rest frames has been applied in refs. [19-22] to 
Lagrangians depending on the acceleration and also on higher derivatives of 
the velocity. In the following we develop a more general formalism which is 
not restricted to rest moving frames. 

The geometry based on the space S of the local Lorentz frames permits 
a formulation of a wide class of (possibly nonlocal) field theories [23] and 
particle theories. In the present article we concentrate our attention on the 
innovative particle models, without a specific consideration of fields differ- 
ent from classical electromagnetism and Einstein's gravitation. The general 
formalism, however, is ready for the introduction of more general classical 
fields. It has been remarked [24, 25] that this geometry permits an elegant 
treatment of the upper bounds to the acceleration and the angular velocity 
of the frames and of the (possibly extended) associated particles. 

An important issue is the description of the motion of a maximally accel- 
erated particle in a Lagrangian or Hamiltonian formalism, since it is necessary 
for the construction of a quantum theory. It has been treated in refs. [21,22], 
by means of the idea, already mentioned above, of associating to the par- 
ticle a moving local Lorentz frame. As it looks natural, rest frames of the 
particle were chosen, but we shall see in the following that this choice is too 
restrictive. 

The report [26] contains a treatment of various kinds of test particles by 
means of Lagrangians that depend on the acceleration and on the angular 
velocity and, for dimensional reasons, must contain a parameter A that can be 
interpreted as a fundamental length. This report also contains a treatment of 
a pre-symplectic formalism, which is more powerful for the treatment of some 
models. In the present article we concentrate our attention on Lagrangians 
not depending on the angular velocity of the moving frame. 

In section II we summarize the relevant aspects of the geometry of the 
fibre bundle S and of its generalizations. In section III we develop a general 
Lagrangian dynamics for a point moving in this space. In section IV we 
introduce the additional assumptions concerning a test particle subject to an 
Einstein-Maxwell field. In section V we use the Noether theorem to find the 
conserved quantities corresponding to various kinds of symmetry properties 
of the system. 

In section VI we give a detailed discussion of the relation between the 
particle and the local reference frames associated to it. We show that, if the 
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Lagrangian depends on the acceleration, the center of mass cannot coincide 
with the origin of the local frame, where the electric charge is assumed to 
be concentrated. In this situation, the concept of acceleration is ambiguous 
and, besides the proper acceleration of the point charge, we define a new 
different quantity, called pseudo-acceleration. 

As a first example, in section VII we give two equivalent treatments of 
the ordinary spinless charged masssive particle. In section VIII we apply the 
general methods developed in the preceding sections to the model treated in 
refs. [21,22], based on rest frames. The model has upper bounds to the veloc- 
ity and to the proper acceleration of the charge, but the energy-momentum 
four- vector may be spacelike. Moreover, even in the absence of external fields, 
the solutions present an unstable behaviour (run-away solutions) when the 
initial conditions are slightly modified. 

In section IX we introduce a different model based on zero-momentum 
frames. It has the correct energy-momentum spectrum and an upper bound, 
which does not concern the acceleration, but the pseudo-acceleration, requir- 
ing a slight reformulation of the maximal acceleration principle. 

In sections X and XI we study further this new model by considering a 
particle in a constant electromagnetic field and we show that the solutions 
have good stability properties. In some circumstances, the velocity of the 
charge may exceed the velocity light, but, since it rotates with a very large 
frequency, this effect is hardly observable. In section XII we consider a 
particle in a simple curved space and, even if the curvature is small, we find 
a serious form of instability. 

In section XIII we analyze another model, with a physical energy-momentum 
spectrum, and the charge moving slower than light. However, the accelera- 
tion and the pseudo-acceleration have no upper bounds. In section XIV we 
draw some general conclusions and we suggest some remarks on the quanti- 
zation. 

II The bundle of Lorentz frames. 

In the absence of an electromagnetic field, a powerful method to treat the 
dynamics of a (possibly extended) particle is to consider, for each value of 
the parameter r, a local Lorentz frame (tetrad) s(r) in such a way that 
a "distinguished point" of the particle (not necessarily the center of mass) 
coincides with the origin of the frame (sometimes called a moving frame). 
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We indicate by S the ten-dimensional bundle of the Lorentz frames of the 
space-time and we describe the motion of the particle by means of the line 
r -> s(t) G S. 

Particle dynamics in this space has been introduced by Kiinzle [12] (see 
also [13]) in the framework of the presymplectic formalism [14]. A treatment 
based on the balance equations for the densities of energy, momentum and 
relativistic angular momentum, strongly influenced by Dixon's work [18], is 
given in refs. [15-17]. In the next section we develop a geometric Lagrangian 
formalism which avoids the introduction of coordinate systems in the space 
S. It is more general than the Lagrangian formalism used in ref. [21], because 
s(t) is not assumed to be a rest frame. 

The elements of S are orthonormal tetrads {eo, . . . , 63} of four-vectors 
in the pseudo-Riemannian space-time manifold M.. We assume that M. is 
time-oriented and eo belongs to the future cone. We consider ten vector 
fields A , . . . ,A 9 in the manifold S. The fields A , . . . , A 3 generate parallel 
displacements of the tetrads along the directions of the tetrad vectors, A 4 = 
A\23], A 5 = ^4 [31], Aq = A\i2] generate rotations around the spatial vectors of 
the tetrad and A 7 = A 10 ], A 8 = Apo], A 9 = A 30 ] generate Lorentz boosts 
along the same spatial vectors. The latin indices i, k, j, I, m, n take the values 
0, . . . , 3 and the greek indices a, (3, 7 take the values 0, . . . , 9. We assume 
A[ik] = A[ki\ an d, when necessary, we use the square brackets to indicate that 
an antisymmetric pair of latin indices stands for a greek index. 

The vector fields A a can also be considered as first order differential 
operators and their commutators (Lie brackets) can be written in the form 

[A a ,A,}=F2,A,. (1) 

The quantities = —Fj a , are called structure coefficients and in the ab- 
sence of gravitation they are the structure constants of the Poincare group. 
In the presence of gravitation, are the components of the curvature 
tensor and F- k are the components of the torsion tensor, which vanishes in 
Einstein's theory. 

The ten dimensional manifold S has a structure of principal fibre bundle 
[27] with base M. and structural group SO" 1 (1,3), but the details of this 
structure are not needed in the next section. The only relevant feature is 
that S is a n-dimensional different iable manifold in which n differentiable 
vector fields A a are defined and that, for each point s G S, the vectors A a (s) 
form a basis of the tangent space T S S. As a consequence, we can identify 
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all the tangent spaces with a single n-dimensional vector space T, which, in 
the absence of gravitation and other external fields, is the Lie algebra of the 
Poincare group. 

In the absence of gravitational field, S is the bundle of the Lorentz frames 
of the Minkowski space-time and the orthochronous Poincare group V acts 
freely and transitively on S. We choose a fixed frame s and represent uni- 
vocally all the other frames in the form s = gs with g e V. In this way we 
may identify S and V . The vector fields A a generate the left translations of 
V. 

We shall also use the differential 1-forms defined by 

^(A a ) = € 

The vectors ou^(s) provide a basis in the cotangent space T*S. 

More details on this kind of geometry can be found in ref. [23]. In this way 
we can also treat a large class of external fields, including torsion fields [1], 
and, for n > 10, gauge fields [28], including electromagnetic fields [29], as it 
is explained in section IV. 



(2) 



III The Lagrange equations. 

We write the tangent vectors of the curve r — > s(r) G S, which describes the 
motion of the frame, in the form 

^ = b a A a (s(r)) e T (3) 
and we consider the action principle 

5 r L{b a , s) dr = 0. (4) 

where the Lagrangian L is an homogeneous function of degree one of the 
"velocities" b a . 

In order to derive the dynamical equations, we consider a family, parametrized 
by e, of varied trajectories (e, r) — > s(e, r) with the usual properties 

s (0, r) = s(r), s(e, n) = s(n), s(e, r 2 ) = s(r 2 ). (5) 

We put 

^ = a"A a . (6) 
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If 4>(s) is a different iable auxiliary function, we have 

^ = ?M = VAJ, (7) 

' >0( ' S) a a A Q (6%0) = 6°A a (a%0). (8) 



dedr 

From the last equality we obtain 



% A a lP)Ap - {b a Ay)A p + a a b p [A a , A p }) = 0, (9) 



namely 

'dW da 1 
de d 



and, since is arbitrary, 



and finally (disregarding higher order terms in e) 

S "° = ' = £ (IF - °^) ( _ = * - a "^)-o' < 12 > 

By means of the last formula, one can write, performing a partial inte- 
gration, 

5 V T l Ldr = e Vrl (fr (¥ - g/3&7 ^ 7 ) + = 

- + A « L ) a ° dT+e [^ a Z ■ (13) 

The last term vanishes as a consequence of the conditions (JHJ) and, considering 
that a a is an arbitrary function of r, we obtain the Euler-Lagrange dynamical 
equations 

p a = ^ = Wp pF P a - A a L, (14) 

where 

P* = " (15) 

If we assume that L does not depend directly on s, the last term in eq. 
()14|) is not present and this equation is exactly the one obtained in ref. [16] in 
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the pole approximation by integrating the balance equations (the quantities 
F2p are defined there with a different sign). 

The momenta p a defined by eq. (fT5j) are homogeneous functions of degree 
zero of the velocities b a , namely they depend only on the ratios 
(and possibly on s). It follows that they must satisfy at least a primary 
constraint [30]. Other primary constraints may exist and we write them in 
the form 

$ p (p a ,s) = 0, p = 0,...,m-l. (16) 

Since L is an homogeneous function, the Euler theorem gives, taking eq. 
(J15|) into account, 

L + Pa b a = (17) 
and by differentiation we obtain 

b a dp a + (A a L) to a = 0. (18) 

The differentials dp a and the forms u a are arbitrary, apart from the con- 
straints 

— M Va + (A Q $ P ) to a = 0, p = 0,...,m-l (19) 
and from eq. (|18p we have 

r = Zcg, (20) 

A a L = J2 aPA «®p- ( 21 ) 
p 

Eq. (J2Uj) inverts, as far as possible, eq. (JT3J). At this level, the func- 
tions a p (r) are arbitrary and, if they are not determined by the dynamical 
equations, they parametrize the gauge transformations of the system. There 
is at least one kind of gauge transformations, namely a redefinition of the 
parameter r. 

If, according to the dynamical equations (JHI), the quantities $ p are not 
conserved, the conditions $ p = determine partially the functions a p or give 
rise to secondary constraints. Some simplification can be obtained from the 
following consequence of the equations derived above 



b a (p a + A a L) = b a b%F^ a = 0. (22) 
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It follows that if a 7^ and the secondary constraints $ p = for p = 
1, . . . , m — 1 are satisfied, $o is conserved. In particular, if there is only one 
primary constraint, it is conserved and there are no secondary constraints. 

If one is able to express all the quantities b a and p a as smooth functions of 
of s and p a (satisfying some constraints), we say that the equations of motion 
are in normal form and one can apply the theorems on the (local) existence 
and uniqueness of the solutions. In fact, if one introduces a local system of 
coordinates q a in the space S, it follows from eq. that the quantities q a 
are smooth functions of the quantities b a . 

IV Einstein- Maxwell fields. 

Now we consider a test particle moving in an Einstein gravitational field and a 
Maxwell electromagnetic field. In order to describe the electromagnetic field, 
we adopt the procedure indicated in refs. [16,28,29], namely we introduce a 
principal fibre bundle S with base M. and structural group SO\l, 3) x U{1) 
which includes the electromagnetic gauge group. Then the manifold S has 
dimension n = 11 and we have to introduce a new vector field, that we 
indicate by A, (in order to avoid a two-digit index a = 10), which generates 
the global electromagnetic gauge transformations. We use the notation 

b a A a = UAi + -b [ik] A [ik] + b*A„ b [ik] = -b [ki] . (23) 

The eleven-dimensional manifold S can also be considered as a principal fibre 
bundle with structural group U(l) and base S , the ten-dimensional bundle 
of the Lorentz frames considered up to now. We shall consider later the 
projection S — > Sq along the fibers generated by the vector field A,. 

The structure coefficients Fuup coincide with the structure constants of 
the Poincare algebra. They can be written in the form 

F mm = ^gkjSf - s^gijSf - 6r 9kl 6j + 6^ 

-5?g kj 5r + SZgyST + 5^g kl 5f - 5 n k gu5f, (24) 
F[ik]j = fykj - 5 l k g lv (25) 

Fmm = o. F m = °- ( 26 ) 

We also assume that the torsion FL vanishes. 
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The structure coefficients F' k = F ik represent the electromagnetic field 
strength and F@ m = F* ik ^ a = 0. The electromagnetic interaction Lagrangian 
is 

Li = eb', (27) 

where e is the electric charge. 

By means of these equations, we can write eq. (|14|) in the more explicit 
form 

Pi = -bi k Pk + F u Fi = b k G lk - AiL, G lk = \p m F [ i ] + eF lk} (28) 

P[ik] = hPk - b k Pi - Wp jk - b^pij - A [ik] L, (29) 
p. = 0, p. = -e, (30) 

where is the Riemann curvature tensor which describes the gravitational 
field. As we discuss in the next section, the quantities p 1 and —pi k are inter- 
preted as the energy, the momentum and the relativistic angular momentum 
of the particle, measured in the local reference frame s. The quantity —p. is 
the conserved electric charge. 

It is often convenient to use a three-dimensional vector formalism. We 
introduce the vectors 

b = (b\b 2 ,b 3 ), b' = {bW,bW,bM), b" = (31) 

P=-(Pl,P2,Ps), P'= -(P[23],P[31],P[12]), P" = -(P[10],P[20],P[30]), (32) 



f = -(b°)-\F 1 ,F 2 ,F 3 ), (33) 

E = (Foi, F 02 , Fos), B = (F 32 , Fi 3 , F 2 i), (34) 

E = (G*oi, Gq2, G03), B = (G32, G13, G21), (35) 
and we obtain, if L does not depend on s, 

p = -b" ■ p + b ■ f , (36) 

p = -b'xp- p b" + 6°f , (37) 

p' = -bxp-b'xp'-b"x p", (38) 

p" = PQ b - b°p - b' x p" + b" x p', (39) 

f = E+(6°) _1 b x B. (40) 
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The dimension of the phase space, namely the number of parameters 
necessary to define the initial conditions, is given by 2n = 22 minus the 
number of primary and secondary constraints (including the constraint (|30p) 
minus the number of the arbitrary gauge parameters (including the usual 
electromagnetic gauge transformations). 

It is instructive to consider, besides the flat Minkowski space-time, the 
spaces of constant curvature 

FH l] = P (6i8l-5i5l). (41) 

M. is a de Sitter spacetime if p > 0, and an anti de Sitter spacetime for 
p < 0. From eq. (|28Jl we obtain 

G ik = -pp[ik] + eF ik , (42) 

or, in the vector notation, 

E = -pp" + eE, B = -pp + eB. (43) 



V Noether's theorem. 

In order to treat the connection between symmetries and conservation laws 
(Noether's theorem), we consider a vector field 

Y(s) = a a (s)A a (s) (44) 

and the corresponding one-parameter diffeomorphism group exp(eY), which 
transforms the trajectory s(r) into the trajectories s(e, r). If this transfor- 
mation does not change the action, the expression (|13|) vanishes and, taking 
the dynamical equations into account, we obtain the conservation law 

o^Va = constant. (45) 

We assume that L does not depend on s and consider two kinds of appli- 
cations of this general theorem, which use two different kinds of symmetry 
properties. In the first case we require 

% Ap] = (a a F2p - Apa?) A, = (46) 
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and from eq. (|12|) we see that 6b a = 0, and therefore 5L = 0. Note that 
the validity of the conservation law (|45|) does not depend on the form of the 
function L(b a ). 

In a second case we assume eqs. (|2~4ll26j) and consider the infinitesimal 
Lorentz transformation generated by Aua. From eq. (|T2*|) we obtain 

5b a = -eF$ kh b\ (47) 

namely an infinitesimal Lorentz transformation of the quantities b a . If the 
Lagrangian is invariant under this transformation, the quantity p^k] is con- 
served. In particular, if the Lagrangian is a Lorentz scalar function of b a all 
the six quantities pw are conserved. If the Lagrangian is only a rotational 
scalar, only the three quantities p[ rs \ with r, s = 1,2,3 are conserved. These 
conservation laws depend on the invariance properties of the Lagrangian and 
on the form of the structure coefficients Fff^- They are also valid in the 
presence of arbitrary gravitational and electromagnetic fields. 

Now we consider the simplest application of the first kind. As we have 
seen in section II, in the absence of gravitational and electromagnetic fields, S 
is the bundle of the Lorentz frames of the Minkowski space-time and choosing 
a fixed frame s, can be identified with the orthochronous Poincare group 
V. The vector fields A a generate the left translations of V, but one can 
also introduce the vector fields A a , which generate the right translations, 
interpreted as Poincare transformations of the fixed frame s. They commute 
with A a and are given by 

Y = A a (g) = DP a (g)Ap(g), (48) 

where D(g) is the adjoint representation of V, which has the property 

A a D%(g) = -F^ s D s y (g). (49) 

It follows from the Noether theorem that the quantities 

p a = D^ a {g{T))p p {r) (50) 

are conserved. Since they are defined starting from the symmetry under 
spacetime translations and infinitesimal Lorentz transformations of s, they 
have to be interpreted as the components of the energy-momentum and the 
relativistic angular momentum measured in the fixed frame s. 
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We indicate an element of V by (A, x) = (A, 0)(1, x), where x is the trans- 
lation four-vector and A is the Lorentz matrix acting on the contravariant 
components of the four-vectors. If we use the explicit form of the adjoint 
representation, we obtain the equations 

pi = A k ip k , P[ik] = A^i^kPiji] - Xip k + x k pi, (51) 

which show that the quantities Pp{r) are correctly interpreted as as the com- 
ponents of the energy-momentum and the relativistic angular momentum 
measured in the moving frame s(r), as we anticipated in the preceding sec- 
tion. The mass /i and the spin a of the particle are given by the familiar 
Poincare invariant expressions and take the same form when written as func- 
tions of ppij) or of Pfiij). In particular we have 

/i 2 = pipt, a 2 /i 2 = S k S k , S k = 2- 1 e k » l p [ij]Pl . (52) 

This interpretation is also valid if an electromagnetic field is present. In 
fact, from eqs. ()31j ) -P0 j) we see that the influence of the electromagnetic field 
on the derivatives p a is just the one expected for a point charge centered 
at the origin, if the quantities p a represent the components of the kinetic 
energy-momentum and relativistic angular momentum. This problem is also 
discussed in ref. [20]. 

The exact definition of the energy-momentum and of the relativistic an- 
gular momentum of an extended particle in a curved spacetime is somehow 
ambiguous [16,18] and in this case the interpretation of p a can be considered 
as a definition. 



VI Frames and particles. 

In order to interpret the solutions of the dynamical equations, we have to 
clarify the connection between the Lorentz frames s(r) and the physical par- 
ticle. We indicate by n : S — > M. the projection which associates to the 
tetrad s G S its origin x = 7r(s) G M.. The projection of the curve r — > s(r) 
is the world line r — > 7t(s(t)) = x(t) in the the spacetime manifold M.. 
From eq. 0, remembering that the vector fields A$ generate the parallel 
displacements along the four-vectors e^r) of the tetrad s(r), we have 

. {r) = d^r) = 6 , (r)ei(r) _ (53) 
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If r is the proper time, namely if bib 1 = 1, the quantities b l are the components 
of the four- velocity of the origin of s(r), with respect to the same frame. 

The tetrad s(r + dr) differs from the parallel transported tetrad by an 
infinitesimal Lorentz transformation with parameters b^dr. This means 
that the covariant derivatives of the tetrad four-vectors are given by 



dr 

and we obtain the formula 



= -6>( T )e fc (r), (54) 



Dt(t) 

a(r) = — ^ = b% - b% k e k = (V + b l k b k ) ei . (55) 

GST 

If t is the proper time, a is the covariant acceleration four-vector and its 
components in the moving frame s(r) are given by 

a* = V + b\b k . (56) 

A rest frame is defined by the condition b = and the four-velocity of 
the origin of the frame is equal to the tetrad vector eo, which is timelike by 
definition. If r is the proper time, we also have b° = 1, and dr is a time 
measured in the rest frame. We have 

a =£p. = bWe r , a° = 0, a = b". (57) 
dr 

The vector a represents the acceleration of the origin of a rest frame measured 
in the same frame. If the parameter r is arbitrary, we can write the more 
general formula 

a=(& )-V. (58) 

We can introduce more specific moving frames by requiring the stronger 
conditions 

b = 0, b^ = b im = 6 [3i] = o. (59) 

In a Lagrangian model they can be imposed by introducing six Lagrangian 
multipliers. In this case s(r) is a Frenet frame and the quantities k\ = b^ m \ 
hi = b^ 12 \ k 3 = &1 23 1 are the geometric invariants of the world line, which 
depend on the acceleration and its first and second derivatives [19,21]. The 
Lagrangian depends on these invariants as in the models considered in ref. 
[19]. 
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Now we consider frames defined by non purely geometric, but dynamical 
requirements. The condition p = defines a zero-momentum frame. It 
implies that the energy-momentum is a timelike four-vector. In the case of 
a free particle, the center of mass is at rest in this frame, but this is not true 
in general. In a zero-momentum frame p' represents by definition the spin of 
the particle. One can still consider the vector (|58|) but it does not represent 
the acceleration of the origin x{r) and not even the acceleration of the center 
of mass. We call it the pseudo-acceleration. From eq. ()37|) we have 

a = p H. (60) 

Note that po in a zero- momentum frame is the invariant mass and this formula 
coincides with the Newton formula, valid for a point particle. If the particle 
is extended, this formula remains valid if we replace the acceleration by the 
pseudo-acceleration. This may be considered as the physical meaning of the 
pseudo-acceleration. 

In order to describe the position of the center of mass, we have to intro- 
duce for each tetrad a system of coordinates in a suitable open set of M., 
for instance a system of normal coordinates. In the absence of gravitation, 
namely in special relativity theory, one can associate to every local frame s a 
Lorentzian coordinate system in the flat Minkowski spacetime and the space 
coordinates of the center of mass at zero time in this frame can be written 
in the form [31] 

y = -PoV- (6i) 

The condition p" = means that the trajectory of the center of mass crosses 
the origin. A frame with this property is called a central frame. This defi- 
nition can be extended to the case in which gravitation is present. From eq. 
(I39|) we see that, if the spin is zero, a zero-momentum central frame is also 
a rest frame. 

The definition p" = is not Lorentz invariant and it is often replaced by 
the condition 

P[ik]P k = 0, (62) 

proposed by Dixon [18]. A frame with this property is called a Dixon frame. 
It can always be obtained from a central zero-momentum frame by means of 
a Lorentz transformation. 

The models treated in the following sections concern extended particles 
which contain a point charge. In this case, instead of working with central 
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or Dixon frames, it is convenient to assume that the charge lies at the origin 
of s(t). More in general, if the charge is not pointlike, one can require that 
the electric dipole moment with repect to the frame s(r) vanishes. 

VII An ordinary spinless particle. 

It is instructive to consider first some Lagrangians that do not contain b' 
and b" and describe a pointlike spinless particle with mass m and charge e. 
The simplest Lagrangian of this kind is 

L = -mb° + eb\ (63) 

Since it is linear in the "velocities" b a , all the momenta p a are fixed by 
eleven primary constraints, namely we have po = m, p. = — e and the other 
momenta vanish. As a consequence, the spin vanishes and the center of mass 
coincides with the charge. From the dynamical equations (|36H39j) we obtain 

b = 0, ma = m(6 ) _1 b" = f = eE. (64) 

We see that the frames s(r) are zero- momentum rest central frames and 
that the acceleration is given by the usual formula. The time evolution of 
the quantities b°, h' and b* is not determined by the dynamical equations, 
the model has gauge invariances described by five parameters and the phase 
space has dimension six. If we impose the gauge fixing condition b' = 0, the 
frame is Fermi- Walker transported and if we assume b° = 1 the parameter r is 
the usual relativistic proper time. If E = 0, the frame is parallel transported 
and the world line of the particle is a geodesic. 

In view of the following developments, it is interesting to treat with more 
detail a charged particle in a constant electromagnetic field in the absence 
of gravitation. The field is invariant under spacetime translations, but it 
changes if a Lorentz transformation is applied to the frame. If b° = 1, we 
have 

E = b"xB-b'xE, B = -b" xE-b'xB. (65) 
Of course, the Lorentz invariants 

I = ||E|| 2 - ||B|| 2 , J = E B (66) 

are constant. 
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The derivatives (|65|) vanish if we choose the gauge fixing condition 



b' = -em -1 B. (67) 

The components of the field and all the quantities b a are constant and we 
say that this is a "stationary" solution. It follws that the projection of 
the trajectory on the ten-dimansional space So, identified with the Poincare 
group V, is given by a one-parameter subroup, namely 

g(T)=exp(rb a A a )g(0), (68) 

where A a form a basis in the Lie algebra of V and exp(.) is the exponential 
mapping of the group [32]. 

We also consider the Lorentz invariant Lagrangian 

L = -m((b ) 2 -\\b\\ 2 ) 1/2 + eb'. (69) 



We obtain 

Po = mb° ((6 ) 2 - ||b|| 2 )~ V2 , p = m ((6 ) 2 - ||b|| 2 )~ 1/2 b, (70) 
p' = p " = 0, p. = -e (71) 

and the constraint 

p 2 -\\p\\ 2 = m 2 . (72) 

It is sufficient to consider the dynamical equations (|37|). since eqs. 
(JSBJ, and (JBTJj) are consequences of the other equations. We see that the 
quantities b°, b', h" and b* are not determined and they are gauge parameters. 
In particular, the Lorentz transformations assume the character of gauge 
transformations. We may fix the gauge partially by assuming that p = 
and we obtain exactly the equations of the preceding model. The dynamics 
of the particle is the same, but in the second model there is more freedom in 
the choice of the moving frame. 

Alternatively, we may require b' = b" = 0, namely that the moving frame 
is parallel transported, and, also assuming b° = 1, eq. (|37|) takes the familiar 
form 

p = f. (73) 
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VIII A model with a maximal acceleration. 



For dimensional reasons, a Lagrangian which is an homogeneous function of 
degree one of the velocities b a and contains a fundamental length A must 
contain, besides b° and b, some of the quantities Ab' or Ab". In the following 
we consider Lagrangians that depend on b", but not on b', and we treat first 
a model with a maximal proper acceleration equivalent to the one treated in 
ref. [21]. 

It is easy to obtain an upper limit for the quantity (6°) _1 b", but, as we 
have seen in section VI, this quantity coincides with the proper acceleration 
only if s(r) are rest frames, namely if b = 0. This constraint is implicit in the 
formalism of ref. [21], but in the framework of the preceding sections we have 
to enforce it by means of the Lagrange multipliers r\. A simple Lagrangian 
of this kind, that does not contain b', is given by 

L = —m ((6 ) 2 - A 2 ||b"|| 2 ) 1/2 + rj ■ b + eb\ (74) 

We obtain the equations 

b = 0, p = r7, p' = 0, p. = -e. (75) 

The first equality means that we are dealing with rest frames and the second 
shows that p can be chosen freely in the initial conditions, even in contrast 
with the usual requirements on the energy-momentum spectrum. The third 
formula does not imply that the spin vanishes, because the spin is defined as 
the angular momentum in a zero-momentum frame, that in general is not a 
rest frame. 

The other momenta are given by 

Po = m&°((6°) 2 -A 2 ||blT 1/2 , 
p " = m A 2 ((b ) 2 - A 2 ||b"|| 2 r 1/2 b" (76) 

and satisfy the scalar primary constraint 

2 \-2|| //|| 2 2 /TTN 

Po - A Hp II = to , (77) 

which, in agreement with eq. (|22|) . is conserved and does not give rise to 
secondary constraints. We consider only solutions with p > m and, as a 
consequence, b° > 0. From the preceding formulas one obtains the relations 

a = (^-V = \- 2 p V = -A- 2 y = A-^HpI 2 + A 2 m 2 )- 1 / 2 p", (78) 
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which show that the acceleration ||a|| has the upper bound A -1 and the 
distance ||y|| of the center of mass from the origin has the upper bound A. 

The dinamical equation (JHUj) follows from eq. ()39|) and the scalar con- 
straint (|77j). Eq. (J38j) is automatically satisfied and from eqs. (pT7|) and (JH9*j) 
we obtain 

p = - Po h" + 6°f -b'xp, p" = -6°p - b' x p". (79) 

There are no secondary constraints. 

As in the simple model considered in the preceding section, the time 
evolution of the quantities b°, W and b* is not determined by the dynamical 
equations and the model has gauge invariances. Since there are five primary 
constraints and five arbitrary gauge variables, the phase space has dimension 
twelve. We can impose the gauge fixing conditions b° = 1, b' = and we get 
the simpler dynamical equations 

p = -A-y + f, p" = -p. (80) 

This system is in normal form. 
We obtain 

^(Po-||p|| 2 ) = -2f'P. (81) 

Even if initially the four-momentum is time-like, by applying a suitable force 
f , it may become space-like. This means that one cannot simply discard as 
unphysical the states with space-like four-momentum. 

We see from eqs. (J2EJ), and PU|) that, when the particle is subject 
to gravitational fields only, if p' = p" = we have f = 0. It follows that 
p" = b" = p = is a solution which describes an ordinary spinless particle 
with mass po = m moving according to the laws of general relativity. We 
show below that, even in the absence of external fields, there are many other 
solutions, as it is expected, since the phase space has dimension larger than 
six. 

It is shown in refs. [21,22] that, in space-times of constant curvature and 
in the absence of other external fields, the dynamical equations can be solved 
exactly. In fact, from eqs. PHjl and (jHUj) we obtain the equation 

p"=(p+A- 2 )p". (82) 

For p > — A~ 2 it has exponential solutions and for p < — A~ 2 it has peri- 
odic solutions. Considering with more detail the case p = 0, namely a flat 
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spacetime, the general solution has the form 

p" = uexp(A" 1 r) + vexp(-A _1 r). (83) 

The acceleration ||a|| and the distance ||y|| of the center of mass from the 
origin of the moving frame are given by eq. (|78|) . For generic initial conditions, 
they approach their maximum values when time increases or decreases, in 
particular 

lim a = A" 1 ||uir 1 u, u ^ 0. (84) 

T— »0O 

For u = v = 0, we obtain p" = a = 0, namely the uniform motion that 
is expected for an ordinary particle in the absence of external forces. This 
solution, however, is unstable, since small changes of the initial conditions 
lead to solutions with a large acceleration in the past or in the future. We 
have 

p 2 - ||pf = m 2 + 4A~ 2 u- v (85) 

and we see that for some solutions the energy-momentum four- vector is space- 
like. The accelerated solutions recall the run-away solutions of the Lorentz- 
Dirac equation [33], with the important difference that in our case we are 
considering test particles and disregarding the radiated energy. 

In a flat spacetime one can easily compute the motion of the origin of 
s(r), namely of the point charge associated with the extended particle. In 
the generic case, it approaches soon the velocity of light with respect to a 
fixed frame. Of course, the center of mass follows a linear world line and its 
distance from the charge becomes large in the fixed frame. In the moving 
frame, however, as an effect of the Lorentz transformation, it remains small. 

If, instead of considering a test particle, we take into account effects 
depending on higher powers of the charge e, it is possible that the radiated 
energy (bremsstrahlung) prevents the unstable behavior described above. In 
the following sections, however, we try a different solution of this problem. 

IX A model with a maximal 
pseudo-acceleration. 

In order to avoid the unphysical energy-momentum spectrum and the unsta- 
ble behaviour of the model described in the preceding section, a reasonable 
attempt is to consider zero-momentum frames instead of rest frames. This 
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can be obtained by dropping the Lagrange multipliers from eq. (J74|) namely 
by assuming 

L = -m((6°) 2 -A 2 ||b"|| 2 ) 1/2 + e6V (86) 
Besides eq. (J3(Jj) we obtain the vector constraints 

p = 0, p' = 0, (87) 

showing that s(r) is a zero- momentum frame and the particle is spinless. 
It is clear that the four-momentum p is necessarily timelike. The equations 
(jTHJ), the scalar constraint (|77j) and eq. (J?HJ) of the preceding section are still 
valid and, in this case too, we consider only solutions with p > m and, as 
a consequence, b° > 0. Eq. (fTSj). however, does not give an upper bound to 
the acceleration, but to the pseudo-acceleration a defined in section VI. 
The relation 

p = m(l-A- 2 ||y|| 2 )~ 1/2 (88) 

shows that p can be interpreted as a confining potential that binds the 
charge to the center of mass. When a force is applied to the charge, ||y|| and 
Po increase and this explains intuitively why the pseudo-acceleration remains 
bounded. 

As in the model of the preceding section, the dynamical equation (jHSl) 
follows from eq. and the scalar constraint (|77j) , eq. (JHS|) is automatically 
satisfied and from eqs. (pT7j) and (J39~j) we obtain 

p h" = b°f, p" = p b - b' X p". (89) 

In this case too we have gauge symmetries. If we impose the gauge fixing 
conditions 6° = 1 and b' = 0, we get the simplified equations 

p" = \ 2 p Q b" = A 2 f = A 2 (E + b x B), (90) 

p" = Pob. (91) 

Note that the four- vector (b°,b), proportional to the four- velocity of the 
charge measured in the frame s(t), is not necessarily timelike and in section 
X we find examples in which it is actually spacelike. It follows that its time 
component in some other fixed reference frame s may vanish or even become 
negative. This means that the charge "moves backwards in time", namely, 
more exactly, the time coordinate with respect to the frame s is a decreasing 
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function of the parameter r. A further discussion will be given in section 
XIV. 

As we have remarked in the preceding section, if there is only a gravita- 
tional field, but no other external fields, the equalities p' = p" = imply 
that f = 0. It follows that the equalities f = p" = b" = b' = b = define 
solutions in which s(r) is a parallel transported rest frame. These solutions 
describe an ordinary spinless particle with mass po = m moving according to 
the laws of general relativity, but in general there are other different solutions 
with p" 7^ 0. In a space of constant curvature p, as we see from eq. fjlHj) . we 
have f = — pp" and if p ^ — A -2 , we have necessarily f = p" = and the 
only solutions are the ones described above. In this case, the phase space 
has dimension six. 

If B ^ 0, from eq. (j9*Uj) we obtain the secondary constraint 

(E - A- 2 p") ■ B = (92) 

We have eight primary constraints, one secondary constraint and five arbi- 
trary gauge variables. It follows that, if there are no tertiary constraints, 
the dimension of the phase space is eight. If B = in the whole space S, 
as in the constant curvature models treated above, eq. (j9*Uj) gives the three 
secondary constraints 

p" = A 2 E (93) 
and the phase space has dimension six, as in the free case. 

X A particle in a constant electromagnetic 
field. 

It is useful to consider, as in section VII, a particle in a flat spacetime with a 
constant electromagnetic field. Some results will be approximately valid for 
slowly varying fields. We consider first the case in which the invariants 
have the properties / > 0, J = and in the frame s(0) we have B = 0. If we 
choose the gauge b' = 0, we see from eq. ()65|) and (J9~U|) that E = B = 0. It 
follows that along the whole trajectory E, p" and b" are constant, B = and 
b = 0. The frame s(r) is a rest frame that moves with a constant acceleration 
given by 

b"= (m 2 + A 2 ||eE|| 2 )- 1 / 2 eE. (94) 



22 



We see that, in this case, the singular manifold defined by B = is invariant, 
namely all the trajectories crossing it are completely contained in it. The 
acceleration remains bounded even if the electric field is very large. As in 
the case of section VII, s(0) determines the initial conditions completely and 
the solution is stationary. 

If in the frame s(0) we have B ^ 0, we must also specify two compo- 
nents of p" in order to describe the initial condition completely. We show, 
however, that there is a finite number of stationary solutions starting from 
s(0) in which all the components of the field and all the quantities b a and p a 
are constant and s(r) is given by the exponential formula (|68|). The other, 
nonstationary, solutions will be treated perturbatively in the next section. 

We assume first that E x B 7^ 0, since the case in which this vector 
vanishes needs a separate treatment. If we require that the derivatives in 
eqs. (|53j) and (JH9*j) vanish, we obtain after some calculation, 



b = aExB, 
b' = p^e (-(1 - a||B|| 2 )B + a(B • E)E) , 
b" = Po x e ((1 - a||B|| 2 )E + a(B • E)B) , 
pl = m 2 + A 2 e 2 ((1 - a||B|| 2 ) 2 ||E|| 2 + a(2 - a||B|| 2 )(E ■ B) 2 ) . (95) 

The second eq. flS9*j) is satisfied if 



f(a) = (1 - a||B|| 2 ) ((1 - a||B|| 2 )(l - a||E|| 2 ) - a 2 (E ■ B) 2 ) - A- 2 e" 2 m 2 a 
= (1 - a||B|| 2 ) (1 - a||B|| 2 - a||E|| 2 + a 2 ||E x B|| 2 ) - A- 2 e - 2 m 2 a = 0. (97) 



This is, in general, a third degree algebraic equation in the variable a, 
that has one or three real solutions. The polynomial function f(a) has the 
following properties 



ap 2 = A 2 e 2 ((l-a||B|| 2 ) 2 + a 2 (B-E) 2 ), 



(96) 



namely, by substituting the expression for p^, 




(98) 
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It follows that f(a) has only one zero in the interval 

< a < min(A 2 e 2 m- 2 , ||E||- 2 , ||B||- 2 ) (99) 
and, possibly, two zeros in the half line 

a > ||B||- 2 (l + A- 2 e- 2 m 2 (||E|| 2 + HBf)- 1 ) . (100) 

If a lies in the interval, ()99|) we have ||b|| 2 < 1, namely the velocity of 
the charge is smaller that the velocity of light. We see below that, for some 
values of the fields, if we choose a solution in the half line (jl00|) . the charge 
may move faster than light. Note that if ||b|| 2 = A 2 ||b' x b"|| 2 > 1, since 
A || b" || < 1, we have A||b'|| > 1, namely the frame s(t) rotates, with all the 
dynamical vector variables, with an angular velocity larger than A -1 . If A is 
very small, one can only observe an average value of the velocity b, that is 
negiglible, since b is orthogonal to the angular velocity b'. 

Note that the fields E and B are measured with respect to the moving 
frame s(r) and if the particle is very fast can be much larger than the fields 
measured in the laboratory frame. Even if we take this remark into account, 
if A is of the order of Planck's length, in all the experimental situations we 
have 

||E||, ||B|| < A _1 e _1 m. (101) 

Under these conditions, the smallest solution of eq. (|97|) is given by the 
approximate formula 

ai ss A 2 e 2 m -2 . (102) 

Note that when A — ^ 0, this solution tends to zero and the other two solutions 
either become complex or tend to infinity. In this limit we obtain, as it is 
expected, the ordinary model of section VII. 

We have to consider the special case in which B ^ and E x B = 0, 
namely E = fcB. We obtain the following conditions for the stationary 
solutions 

b = 0, b' = /i||B|| -1 B, p b" = eE, p 2 = m 2 + A 2 ||eE|| 2 , (103) 

where h is an arbitrary constant that describes just a choice of the rotational 
gauge. Note that eq. (}9"4*j) is also valid in this more general situation. 

It is interestig to see if these equations can be obtained as a limit of eqs. 
(J95J) for || E x B || — > 0. The equation (JDTj) tends to the second degree equation 

(1 - a||B|| 2 )(l - a||E|| 2 - a||B|| 2 ) - A" 2 e~ 2 m 2 a = (104) 
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and one of the three solutions of eq. (J97|) tends to infinity, more exactly we 
have 

a 3 « (||E|| 2 + ||B|| 2 )||E x B||" 2 . (105) 

It follows that for this solution we have ||b|| — > oo, showing that there are 
stationary solutions with the charge moving faster than light. The other two 
(necessarily positive) solutions are given by 

oi, 2 = (2||B|| 2 (||E|| 2 + HBf^fllEH 2 + 2 1| B || 2 + A~V 2 m 2 ± A 1/2 ). (106) 

where 

A = (||E|| 2 + A-V 2 m 2 ) 2 + 4A- 2 e- 2 m 2 ||B|| 2 > 0. (107) 
The corresponding values of the parameter h are 

/» 1 ^=Pb 1 c||B||(a ± (||E|| 2 + ||B|| 2 )-l) 
= e(2p ||B||)- 1 (||E|| 2 + A- 2 e- 2 m 2 ± A 1/2 ). (108) 

Only for these values of the gauge parameter h the solutions ()103|) can be 
obtained as limits of the solutions (J95j) . 

XI Perturbations and stability. 

The next step is to study perturbatively trajectories slightly different from 
the stationary ones described above. Since the infinitesimal perturbations 
satisfy linear differential equations with constant coefficients, we look for 
complex exponential solutions, namely we write 

b a -> b a + 3? (5b a exp(zr)) , (109) 

where the constant quantities z and Sb a are complex. A similar notation 
is used for the other variables. The derivatives with respect to r can be 
replaced by the factor z. We are not interested in perturbations with z — 0, 
which give other already known stationary solutions. For z ^ we obtain 
new nonstationary solutions, but also stationary solutions, modified by an 
infinitesimal time-dependent gauge transformation. 
From eq. (|65j) we obtain 

z5E = b" x 5B - B x 5b" -b'xffi + Ex 5b', 
z5B = -b" x5E + Ex 5b" -b'x<5B + Bx 5b'. (110) 
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Since the invariats (|66|) are not affected by the perturbation, we also require 
the relations 

E-<5E = B-5B, E 5B = -B <5E, (111) 

which, however, are not independent from eqs. (jllUj) . 
From the other dynamical equations we have 

5p" = X 2 e{5E + b x SB - B x 5b), 
z5p" = p Q 5b + b5p - b' x 5p" + p" x 5b', 
5p" = X 2 b"5p + X 2 p 5b", 
5p = X 2 p 3 m- 2 b" ■ 5b". (112) 

After the elimination of the variables 5po, 5p" and 5b, we obtain a system of 
nine homogeneous linear equations in the twelve unknown variables 5E, 5B, 
5b" and 5b'. 

For all the values of z this system has three linearly independent solutions 
which represent rotational gauge transformations. They are given by 

<5E = Exr, 5B = Bxr, 5b" = b" x r, 5b' = b'xr + zr, (113) 

where r is an arbitrary infinitesimal complex vector. Note that unless 

^ = -||b'|| 2 , (114) 

the vector 5b' can take arbitrary values and when we look for physically 
relevant perturbations we can choose the gauge 5b' = 0. We obtain in this 
way an homogeneous linear system of nine equation in nine unknowns, which 
has nonvanishing solutions only if its determinant det M(z) vanishes. If z 
satisfies eq. (I114j) . there is a pure gauge perturbation with 5b' = 0, and this 
implies that det M(±z||b'||) = 0. If we disregard these uninteresting solution, 
the other nonvanishing solutions of the algebraic equation det M(z) = 
correspond to physically relevant perturbations. 

If E x B ^ 0, it is convenient to use the projections of all the vectors 
on the basis formed by the vectors E, B and E x B. Eq. (|111|) permits the 
elimination of two unknowns and we have to calculate the determinant of a 
7x7 matrix. It has almost thousand terms and one has to use a computer 
algebra program [34]. We find that the determinant detM(^) contains only 
the powers z 6 , z 4 and z 2 , and, since we know that it vanshes for z 2 = — ||b'|| 2 , 
we can find the other, physically relevant, solution of the eq. detM(z) = 
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by means of Ruffini's rule. With an appropriate choice of the variables, the 
computer calculation gives the unexpectedly simple result 



-2 , .2 \-4 2 



Z = UJ 



A- 4 pol|eB|p + AX~ 2 m 2 p~ 2 > 0, (115) 



where p§ i s given by eq. (J95J) . We see that z = ±iu is pure imaginary 
and this means that the stationary solutions are stable with respect to linear 
perturbations. The perturbations have a frequency u/2ir, that when ||eB|| — > 
0, tends to infinity explaining why a degree of freedom disappears for ||eB|| = 
0. If we use eq. (J101J) . we see that the frequency u/2tt is much larger than 
the Planck frequency A -1 . 

In order to extend our results to the case in which E x B = a separate 
discussion, that does not require computer algebra, is needed. From eqs. 
(PUD (fTTTjl and (HI2J) we find 

B-5E = B-5B = B-5p" = B-5b = 0, 5p = 0. (116) 

As in the general case, we choose 5h' = 0. 

After some calculations we obtain the relation 

z{5B - \ 2 kepQ 1 \\B\\ 2 5v") = -ti x (SB - X 2 ke P Q 1 \\B\\ 2 5p") (117) 

and if z 2 ^ — ||b'|| 2 we have 

= A" 2 (/fce)- 1 pol|B||~ 2 5B. (118) 

By substituting this formula and eq. (|1U3|) into eq. (JTTUJ), we obtain a 
linear homogeneous system of the form 

^b)=P^(bxS). , 119 , 

where the matrix A is given by 

A ~ \ -PoleEH -/i + A-^olleBU- 1 J' 

Since the vectors 5E and 5B must be perpendicular to B, by iterating this 
formula we obtain 

*(S )--*(£)■ <-) 
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One can easily show that this linear system has nonvanishing solutions if 

z 2 = -(/ii, 2 - h) 2 , (122) 

where the quantities hip are given by eq. (|1U8|) . From the same equation and 
eq. ()115j) we obtain the identity 

h 2 -h 1 =u. (123) 

We see that, in this case too, the stationary trajectories are stable, but the 
frequency of the perturbations depends on the gauge parameter h. There are 
two independent perturbations with different angular velocities and only their 
difference is invariant under the rotational gauge and physically relevant. If 
we choose h = hi i2 , in both cases we obtain the solution z 2 = and the 
solutions given by eq. (|115|) . namely the limits of the solutions found for 
ExB^O. 



XII A particle in a simple gravitational field. 

In order to understand the properties of the model in a gravitational field, 
it is instructive to consider, besides the spacetimes with constant curvature 
we have already examined in section IX, another simple model, namely the 
stationary Einstein cosmological model [35] with spatial curvature p. The 
curvature tensor in the frame s(r) takes the form 

= p(6{6i - 8{5\ - b\v k v l + 6{v t v l + 5\v k v> - 8{v^), 

v k v k = 1, p > 0, (124) 

where, in a cosmological interpretation, v — (vq, v) is the four- velocity with 
respect to the frame s G S of a privileged frame in which the matter and 
radiation distribution is isotropic. 
From eq. (|28|) we obtain 

G ik = -p(p[ik] + v 3 p [:j i]V k - v^p^^Vi), (125) 

or, in the vector notation, 

E = p(||v||V - (v ■ p», B = pv v x p". (126) 
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Eqs. (JUD (JHD give the result 

{\- 2 p- 1 - ||v|| 2 )p" + (v ■ p> = po^op" x (v x p"), (127) 

and the constraint (|9*2*|) is identically satisfied. We assume that ||v|| 2 < 
A~ 2 p _1 . If v x p" = we must have p" = 0, otherwise we obtain 

v • p" = -p,v \\- 2 p- 1 - ||v|| 2 ), p = A- V(v ■ p"). (128) 

The component of p" normal to the vectors v and p" is not determined by 
the dynamical equations, namely there is a gauge symmetry. We fix the 
gauge by means of the condition 

p" . v x p" = 0. (129) 

One can show that the covariant spacetime derivatives AiV k vanish and 
it follows that the derivatives of v k with respect to the parameter r are given 
only by the boost of the Lorentz frame, namely 

v = -v b" = -\- 2 p^v p", v = -(v ■ b") = -\- 2 p l (v ■ p"). (130) 

From the preceding equations it follows 

< ||v x p"|| 2 = \\pl - m 2 )||v|| 2 - A\ 2 (p ) 2 , (131) 
dw 

— = 0, w = p v > m. (132) 
dr 

The conserved quantity w is the energy measured in a privileged frame. We 
also obtain 

^ logpo = -A" V^A-V 1 + 1 - w 2 m 2 p^). (133) 

If w is not extremely large, the right hand side is negative, Po{t) has a 
maximum for r = f and violates the inequality Po > m outside a given 
interval. A stronger limitation is given by eq. ()131|) . It is clear that r cannot 
vary in the whole real line. 

For a detailed discussion, it is convenient to introduce some simplifica- 
tions. We disregard higher powers of the adimensional quantity A 2 po, which 
is extremely small in all the interesting cases and we also disregard higher 
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powers of r — r. We put Po( r ) = Po, v o( T ) — v o, v ( r ) = v an d we obtain the 
approximate solution 

p (r) » p (l - 2~ X A- W(t - f) 2 ), (134) 

Po(r) « -p A~ W(r " r). (135) 
From eq. ()131|) we see that r must belong to the interval 

|r - f | < 2~ 1 At « (1 - m 2 po 2 ) 1/2 AVo||v||. (136) 

This inequality shows that it was consistent to disregard higher powers of 
t — t together with higher powers of X 2 p. In this interval we have 

< Po ~ Po(r) < 2- 1 A 2 pp (l - m% 2 )||v|| 2 (137) 

and the quantities p , \\p"\\, v , and v can be considered as constant. 

It follows from eq. f)131j) that at the end points of the interval ()136|) p" 
is parallel to v and by taking eq. (|129|) into account, we see that the vector 
p" moves approximately on an half circumference in a plane containig the 
vector v. From eq. ()128|) we see that the projection of p" in the direction of 
v is constant and it follows that the component of p" normal to v tends to 
infinity at the end points. 

The connection between r and the "cosmic" time t, measured in the 
privileged rest frames defined by v = 0, is 

^ = Vl U = Vo -vb= (A-V 1 + IK" 1 (138) 
ar 

and it follows that the interval ()136|) measured in terms of the time t is 

At w 2A(1 -m 2 p 2 ) 1/2 ||v||, (139) 

namely it is, in general, very small. Note that it does not depend on the 
curvature p, but only on the four- vector v that defines the privileged frames. 

The singularities at the end points of the interval ()13ti|) can be avoided 
by using, instead of the parameter r defined by the condition b° = 1, another 
parameter and we can also try to extend the solution outside the interval. 
This is possible only if, generalizing our Lagrangian formalism, we accept 
negative values of b°, namely negative values of the square root in eq. (|75|) . 
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We do not try to give in this paper a physical meaning to this extension, for 
instance in terms of antiparticles. 

In any case, we have shown that, in the presence of an arbitrarily small 
gravitational field, the geodisic world lines, corresponding to the initial con- 
dition p" = 0, are badly unstable, since an arbitrarily small change of p" 
leads to completely different, physicaly unacceptable, solutions. 

XIII Another model. 

In the preceding sections we have examined two models that imply an upper 
bound to acceleration or to pseudoacceleration. However, these models have 
unwanted properties, in particular the energy-momentum or the four-velocity 
of the charge may become spacelike and unstable solutions may appear. In 
order to show that it is not easy to avoid all these problems, we consider a 
model containing a fundamental length, in which both the energy-momentum 
and the four-velocity are necessarily timelike, as in the ordinary model of 
section VII. A simple and natural Lagrangian with these properties is 

L — —m ((b ) 2 - ||b|| 2 - A 2 ||b"|| 2 ) 1/2 + eb\ (140) 
The canonical momenta are given by 

p' = 0, p. = -e, (141) 

p = mb° ((b ) 2 - ||b|| 2 -A 2 ||b"|| 2 r 1/2 , 
p = m((fo ) 2 - ||b|| 2 -A 2 ||b"|| 2 )" 1/2 b 
p" = m A 2 ((b ) 2 - \\b\\ 2 - A 2 ||b"|| 2 )" 1/2 b" (142) 

and satisfy the scalar primary constraint 

(po) 2 -||p|| 2 -A-VT = m 2 , (143) 

which assures that the energy-momentum is timelike. 

We adopt the the gauge fixing conditions b° = 1 and b' = and we obtain 
the formulas 

b = p-tp = (A- 2 ||p"|| 2 + ||p|| 2 + m 2 )- 1 ' 2 ?, (144) 
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b" = A-VV = A _1 (IIp"I| 2 + A 2 ||pf + X 2 m 2 )- 1/2 P" , (145) 

which assure that ||b|| < 1 and ||b"|| < A -1 . Note, however, that s(r) is nei- 
ther a rest frame nor a zero-momentum frame and b" cannot be interpreted 
as the acceleration or the pseudoacceleration of the particle. 

The dinamical equation (j3T?j) follows from eq. ()39|) and the scalar con- 
straint (|143|) . Eq. (J3SJ) is automatically satisfied and from eqs. (I37J) and 
we obtain, 

p = - Po b" + f = -A- 2 p" + f , p" = p h - p = 0. (146) 

We see that p" is conserved. If it vanishes, the dynamical equations coincide 
exactly with the equations of the ordinary model of section VII in its last 
version and there is no upper bound to the acceleration or to the pseudoac- 
celeration. 

XIV Conclusions. 

We have developed a general Lagrangian formalism for the description of test 
particles in gravitational and electromagnetic fields. It is based on a moving 
Lorentz frame associated to the particle and it is useful for the introduction 
of Lagrangians that depend on a fundamental length A and to discuss a 
possible upper bound to the proper acceleration. Other interesting effects 
may be suggested by a detailed discussion of suitable models. 

The test particles introduced in this way present some features charac- 
teristic of extended objects, in particular, additional degrees of freedom are 
present, the charge, assumed to be concentrated at the origin of the moving 
frame, does not necessarily coincide with the center of mass, the energy- 
momentum four-vector is not necessarily parallel to the four-velocity and the 
definition of the acceleration is somehow ambiguous. 

We have considered with more detail three Lagrangians depending on 
the acceleration of the frame, but not on its angular velocity, discussing, in 
particular, the following requirements: 

• The energy-momentum four-vector lies in the future cone; 

• The velocity of the point charge is smaller than the velocity of light, 
namely the four-velocity is timelike; 

• The solutions have good stability properties; 
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• There is an upper bound to the acceleration or to some related quantity 
as the pseudo-acceleration defined in section VI. 

We have not been able to find a model that satisfies all these requirements 
and it is reasonable to conjecture that such a model does not exist within the 
formalism we have considered. Perhaps one should introduce new degrees of 
freedom that have not a geometric interpretation in terms of a moving frame. 

The model deascribed in section VIII, based on rest frames, has an un- 
physical energy-momentum spectrum and has unstable run-away solutions 
even in the absence of external fields. 

The model introduced in section IX, based on the zero-momentum frames, 
has unstable solutions when a nonconstant curvature is present and, under 
some circumstances, the velocity of the charge may become arbitrarily large. 
Nevertheless, it has some interesting properties. It always has a timelike 
energy-momentum four-vector and it has good stability properties if only a 
slowly varying electromagnetic field is present. 

An electric charge faster than light may be embarrassing, but does not 
seem to contradict the general principles of Maxwell's theory. The charge 
describes a singular distribution inside the particle of the electric current 
four-vector, that is not required to be time-like. Moreover, in the absence of 
curvature, the velocity vector rotates very fast and one can only observe its 
average value. 

In a constant electromagnetic field the second model has "stationary" 
solutions, similar to the trajectories of an ordinary point particle with some 
small corrections, but, since the model has additional degrees of freedom, 
more complicated solutions exist, obtained by adding to a stationary solution 
a periodic perturbation with a frequency u)/2n much larger than A -1 . 

In a quantized version of the model, according to Bohr's correspondence 
principle, this frequency is approximately interpreted as the frequency of the 
radiation emitted when an excited state of the particle decays. Since the 
frequency is defined in the zero momentum frame, the energy hu has to be 
interpreted as the additional mass of the excited state. This mass is very large 
and in ordinary conditions there is not enough energy to excite the periodic 
degrees of freedom and they are "frozen", as it happens, for instance, to 
the nuclear degrees of freedom in low-energy atomic physics. These features 
of the quantized model may provide some very simplified insight into the 
connection between ordinary and "transplanckian" physics. 

The model presented in section XIII has no problem with the properties 
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of the energy-momentum and of the four-velocity, but the presence of the 
fundamental length A in the Lagrangian does not impose any limitation to 
the acceleration. 

Of course, a more precise physical discussion requires a formal quanti- 
zation procedure of the models, based on Dirac's treatment of constrained 
systems. In particular, the ground state energy corresponding to the addi- 
tional (internal) degrees of freedom, which cannot be evaluated by means of 
the correspondence principle, could give an unacceptably large contribution 
to the mass. The treatment of many particle states should be given in terms 
of free quantum fields. For the model of section VIII these problems are 
treated in ref. [36] and we shall discuss some other models elsewhere. 
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